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Abstract 



The relativistic nuclear recoil corrections to the energy of the 2p3 

2 

state of hydrogen-like and the (ls)^2p3 state of high Z lithium-like 
atoms in all orders in aZ are calculated. The calculations are car- 
ried out using the B-spline method for the Dirac equation. For low 
Z the results of the calculation are in good agreement with the aZ 
-expansion results. It is found that the total nuclear recoil contribu- 
tion to the energy of the (ls)^2p3 — (ls)^2s transition in lithium-like 
uranium constitutes —O.OQeV. 
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1 Introduction 



Till recently accurate QED calculations of the nuclear recoil corrections to 

atomic energy levels were of interest mainly in connection with high precision 
measurements of the Lamb shift in hydrogen [1,2]. In this case the nuclear 
recoil corrections may be calculated in lowest orders in aZ [a is the fine 
structure constant, Z is the nuclear charge). However recent achievements 
in experimental investigations of highly charged ions [3] require calculations 
of the nuclear recoil corrections right up Z = 92. In the last case the param- 
eter aZ can no longer be considered small and, so, the calculations without 
expansion in aZ are required. 

In our previous paper [4] we calculated the nuclear recoil corrections in 
all orders in aZ for the Is, 2s, 2pi states of hydrogen-like atoms and the 
(ls)^2s, (ls)^2pi states of high Z lithium-like atoms. We found that the 
nuclear recoil contribution to the energy of the (ls)^2pi - (ls)^2s transi- 
tion in lithium-like uranium constitutes —0.07eV and, so, is comparable 
with the uncertainty of the experimental value of the transition energy: 
280.59(10) eV [3]. We found also the nuclear recoil contribution, additional 
to Salpeter's one, to the Lamb shift (n = 2) of hydrogen to be —1.32(6) kHz. 
This result is in good agreement with recent analytical calculations of the 
^{aZ)^ In (aZ) and ^(aZ)^ corrections [5-8], according to which the terms 
of order '^{aZ)^\n{aZ) cancel each other [5,6], the contribution of order 
^{aZf is -0.77 kHz for the 2s state [7] and 0.58 kHz for the 2pi state [8]. 

In the present paper we extend the results of [4] to the case of the 2»3 

2 

state of hydrogen-like and high Z lithium-like atoms. 
The relativistic units h — c— 1 are used in the paper. 

2 Hydrogen- like atoms 

The complete aZ-dependence expressions for the nuclear recoil corrections 
to the energy levels of hydrogen-like atoms were first derived by one of the 
authors of the present paper [9] (a part of the expressions was found ear- 
her by Braun [10]). These expressions consist of three contributions: the 
Coulomb contribution, the one-transversc-photon contribution, and the two- 
transverse-photon contribution. For a state a the Coulomb contribution is 
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given by 

AE, = A£;«+A£;f , (1) 



= (al^la), (2) 

A^f) = ^y_^rfu;5^(a;)(a|[p,V;]G(a; + £„)[p,K]|a), (3) 

where |a) is the unperturbed state of the Dirac electron in the Coulomb field 
of the nucleus, K = — ^ is the Coulomb potential of the nucleus, p is the 
momentum operator, 5+ (a;) = ^('^ + ^O)"-*^, G(oi)) — (oj — H(l — iO))~^ is 
the relativistic Coulomb Green function, H — ap + Pm + Vc . The one- 
transverse-photon contribution is 

, (4) 
D(0)p + pD(0))|a), (5) 



A-E'tr-(l) 










M J- 



dw 5+{uj){a\ ([p, Vc]G{uj + £„)D(a;) 
-DHG(a; + £„)[p,K])|a), (6) 



where 



Dm{uj) = -AiraZaiDimiuj) , (7) 



ai {I = 1,2,3) are the Dirac matrices, Dim{uj) is the transverse part of the 
photon propagator in the Coulomb gauge. In the coordinate representation 
it is 

DA-.r) = -±{^^^^6^ + v.v, ''^P''lf'-^) } . (8) 

47r r cu^r ' 

The two-transverse-photon contribution is 

^Etri2) = ^ j_^duj{a\-D{uj)G{u + ea)T>{uj)\a) . (9) 

The terms A£^^^) and AE^^-^^^ are leading at low Z. These terms can easily 
be calculated by using the virial relations for the Dirac equation [11-13]. Such 
a calculation gives [9] 
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A£;W ^ A£;W + A<) ) = . (10) 

The aZ-expansion of this term is given in [4]. 

The terms A^^^)^ ^E^i^i). and ^Etr{2) are given in the form that allows 
one to use the relativistic Coulomb Green function for their calculation. (As 
was found by Pachucki and Grotch [7], this form is convenient for the aZ- 
expansion calculations as well.) In [4] we transformed these equations to 
the form that is more convenient for calculations using the finite basis set 
methods [14] and calculated them for the Is, 2s, and 2pi states by using the 
B-spline method for the Dirac equation [15,16]. The calculation for the 2^3 
state considered here was carried out in the same way. 

Table 1 shows the results of the numerical calculation for the 2»3 state 

2 

expressed in terms of the function P[aZ) defined by 



AE(^) = A^f + ^E%, + AE,.(.) = (11) 
The functions Pc? Ptr{i), and Ptr{2) correspond to the contributions AE^^\ 

(2) 

^-^tV(i)' ^Etr{2), respectively. For comparison, the Salpeter's contribu- 
tion [17] is 

Ps{aZ) = -0.030017- ^ = -0.30884. (12) 
3 18 

The uncertainties given in the table correspond only to errors of the numerical 
calculation. In addition, there is an uncertainty due to deviation from the 
point single particle model of the nucleus used here. 

To make a more detailed comparison with the aZ-expansion calculations 
we represent the function P{aZ) in the form 

P{aZ) — ai + a2aZ + a3{aZy + a4{aZ)^ (13) 

The coefficients can be calculated from our numerical results for the 
P(a;Z)-function. Such a calculation using the values of the P(aZ)-function 
for Z =1,2,3,5,10,15 gives 

ai = -0.30883, as = 1.040. (14) 
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The coefficient ai is in good agreement with the Salpeter's resuh given by 
the equation (12). The coefficient a2 coincides, within errors of the numer- 
ical procedure, with the corresponding coefficient (02 = 1.047) obtained by 
Golosov et.al. [8] ( the related coefficient from the contribution (10) is equal 
to zero). 

According to our numerical calculations the contribution of the difference 
between AE^'^'' and Salpeter's correction for the state of hydrogen con- 
stitutes 0.42(2) kHz. Analytical result for the ^{aZ)^ correction found in 
[8] gives 0.423 kHz. ( In addition to this correction, in [8] the correction of 
order '^o?{aZY is calculated to be 0.014 2;.) 

Let us consider the nuclear recoil corrections for hydrogen-like uranium. 
According to the formula (10) the first correction is 

AE-g^g = 0.0664 ey. (15) 

The second correction defined by (11) is 

AEg^^ = 0.0014 ey. (16) 

In the next section we use these results to find the total nuclear recoil con- 
tribution to the energy of the 2'pz — 2s transition in lithium-like uranium. 



3 High Z lithium-like atoms 

The complete aZ-dependence expressions for the nuclear recoil corrections to 
the energy levels of high Z few-electron atoms were derived in [18]. These cor- 
rections are the sum of one- and two-electron contributions. The one-electron 
corrections are obtained by summing all the one-electron contributions over 
all the one-electron states that are occupied. The two-electron corrections 
for the case considered here are 

A£;f*^ - i: (a|p|n)(n|p|a), (17) 

A<5 = ]^ E {{a\p\n){n\B{e^-e^)\a) (18) 
+ (a|D(£„ - £„)|n)(n|p|a)} , 
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A^SpJ = E HD(£a-£n)|n)HD(£„-£j|a). (19) 

£n=£ls 

The table 2 shows the results of the calculation of the corrections (17), (18), 
and (19) for the (ls)^2p3 state expressed in terms of the function Q{aZ) de- 
fined by 

AE^nt ^ ^Ei^nt) + ^^M) + ^^M) ^ _^^(«z)2Q(aZ) . (20) 

Here we have taken into account the known non-relativistic limit of this 
correction [19]. Within the ^(ctZ)^ approximation the function Q{aZ) that 
we denote by Q^laZ) is [20] 

QUaZ) = l + + (21) 

For comparison, this function is given in the tabic as well. The functions 
QciaZ), Qtr{i){aZ), and Qtr{2){oiZ) correspond to the corrections Ai?^*"*\ 

, and AEj.^^^) ) respectively In leading orders in aZ they are 

17 1 97 

QciaZ) = ^ + i^Z)'{- + -ln-), (22) 

QMi)(a^) = liaZ)\ (23) 
25 

Q,,(2)(aZ) = -^(«^)'- (24) 

For low Z, in addition to the corrections considered here, the Coulomb 
electron-electron interaction corrections to the non-relativistic nuclear recoil 
contribution must be calculated separately. The main contribution from 
these corrections is of order ^iaZ)^^. 

Let us find the total nuclear recoil contribution to the energy of the 
(ls)^2p3 — (ls)^2s transition in lithium-like uranium. According to our cal- 
culation the term AE^"-^ contributes — 0.0345 el^. Adding to this value the 
one-electron contribution defined by (15) and (16) and using the related val- 
ues for the 2s state from [4] we find 

A^(i,)22p3 - A^(i,)22, = -0.094 

This value is significant enough to be included in accurate QED calculations 
of the transition energy. 
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Tabic 1: The results of the numerical calculation of the one-electron nuclear 
recoil corrections to the 2pa state energy expressed in terms of the function 
P{aZ) defined by equation (11). 



z 


Pc(aZ) 


Pfr(u(aZ) 


Pfr(2)iOiZ) 


P{aZ) 


1 


0.0000 


-0.1385(2) 


-0.1630(3) 


-0.3013(4) 


5 


-0.0001 


-0.1230 


-0.1493(1) 


-0.2724(1) 


10 


-0.0004 


-0.1037 


-0.1338 


-0.2379 


15 


-0.0008 


-0.0841 


-0.1198 


-0.2047 


20 


-0.0012 


-0.0644 


-0.1069 


-0.1726 


25 


-0.0017 


-0.0446 


-0.0950 


-0.1413 


30 


-0.0021 


-0.0245 


-0.0841 


-0.1107 


35 


-0.0026 


-0.0042 


-0.0741 


-0.0081 


40 


-0.0031 


0.0163 


-0.0649 


-0.0517 


45 


-0.0035 


0.0372 


-0.0567 


-0.0230 


50 


-0.0040 


0.0585 


-0.0496 


0.0050 


55 


-0.0044 


0.0802 


-0.0432 


0.0326 


60 


-0.0049 


0.1025 


-0.0379 


0.0597 


65 


-0.0054 


0.1254 


-0.0337 


0.0863 


70 


-0.0058 


0.1490 


-0.0307 


0.1125 


75 


-0.0063 


0.1735 


-0.0289 


0.1384 


80 


-0.0068 


0.1990 


-0.0285 


0.1638 


85 


-0.0072 


0.2256 


-0.0294 


0.1889 


90 


-0.0078 


0.2534 


-0.0319 


0.2138 


92 


-0.0080 


0.2649 


-0.0334 


0.2236 


95 


-0.0083 


0.2827 


-0.0361 


0.2383 


100 


-0.0088 


0.3135 


-0.0422 


0.2625 
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Table 2: The results of the numerical calculation of the nuclear recoil correc- 
tion 

^^(int) ^Yie {lsY2ps state of lithium-like ions expressed in terms of 
the function Q{aZ) defined by equation (20). Q]:^{aZ) is the leading contri- 
bution defined by equation (21). 



z 


QJoiZ) 


Qtr(n(ciZ) 


QtrCl^iciZ) 


Q{aZ) 


QL(aZ) 


5 


1.00036 


-0.00083 


0.00000 


0.99953 


0.99953 


10 


1.00143 


-0.00333 


0.00000 


0.99810 


0.99811 


15 


1.00323 


-0.00748 


0.00001 


0.99573 


0.99574 


20 


1.00574 


-0.01330 


0.00004 


0.99239 


0.99242 


25 


1.00896 


-0.02077 


0.00010 


0.98809 


0.98816 


30 


1.01291 


-0.02989 


0.00021 


0.98281 


0.98295 


35 


1.01757 


-0.04065 


0.00039 


0.97653 


0.97679 


40 


1.02294 


-0.05303 


0.00065 


0.96926 


0.96969 


45 


1.02902 


-0.06704 


0.00101 


0.96097 


0.96163 


50 


1.03579 


-0.08264 


0.00150 


0.95165 


0.95264 


55 


1.01323 


-0.09981 


0.00212 


0.91129 


0.91269 


60 


1.05132 


-0.11854 


0.00289 


0.92988 


0.93179 


65 


1.06001 


-0.13877 


0.00382 


0.91742 


0.91995 


70 


1.06925 


-0.16047 


0.00488 


0.90390 


0.90716 


75 


1.07894 


-0.18356 


0.00608 


0.88930 


0.89343 


80 


1.08897 


-0.20796 


0.00738 


0.87362 


0.87875 


85 


1.09914 


-0.23356 


0.00873 


0.85686 


0.86312 


90 


1.10921 


-0.26019 


0.01006 


0.83896 


0.84654 


92 


1.11313 


-0.27109 


0.01057 


0.83148 


0.83964 


95 


1.11880 


-0.28765 


0.01127 


0.81988 


0.82901 


100 


1.12736 


-0.31561 


0.01224 


0.79951 


0.81054 
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